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ABSTRACT 


The  most  vital  link  in  a  single  cummodity  How  network  is  that  are  whose  re¬ 
moval  results  in  the  greatest  reduction  in  the  value  of  the  maximal  flow  in  the  network 
between  a  source  node  and  a  sink  node.  This  paper  develops  an  iterative  labeling 
algorithm  to  determine  the  most  vital  link  in  the  network.  A  necessary  condition  for  an 
arc  to  be  the  most  vital  link  is  established  and  is  employed  to  decrease  the  number  of 
arcs  that  must  be  considered. 
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SEiCTION  1 
INTRODUCTION 


The  problem  of  removing  arcs  and  nodes  from  a  network  is  a  part  of  network 
theory  that  has  many  Mnportant  and  useful  applications.  One  application  would  be  u 
conflict  situation  where  there  is  a  logistics  or  communications  network  under  attack.  A 
defender  or  user  of  the  system  must  know  which  arcs  arc  most  vital  to  him  so  that  he 
can  reinforce  them  against  attack;  white  the  attacker,  naturally,  wants  to  destroy  those 
arcs  whose  destruction  would  most  affect  the  efficiency  of  the  system.  Another  applica¬ 
tion  would  be  in  helping  the  managers  of  a  highway  system  or  a  transportation  network 
to  determine  the  effect  of  closing  various  links  for  repair,  etc. 

The  problem  addressed  by  this  paper  is  concerned  with  finding  the  most  vital 
link  in  a  single  commodity  flow  network  (N;A1  (directed,  undirected,  or  mixed).  An  arc 
(x,y)  6  A  is  declared  to  be  the  most  vital  link  if  its  value  v(x,y)  is  at  least  as  large  as  the 
value  of  every  other  arc  in  the  network.  The  value  of  arc  (x.y)  is  defined  as  the  difference 
in  the  maximal  flow  values  in  networks  (N;A)  and  j^N;|A-(x,y)j  between  some  source 
node  and  some  sink  node.  Thus,  v(x,y)  reflects  the  reduction  in  t  le  maximal  flow  attain¬ 
able  if  arc  (x.y)  is  removed  from  the  network. 

Wollmer  has  developed  an  algorithm  for  determining  the  most  vital  network  link.^  *  ^ 
This  algorithm  has  been  employed  by  Durbin  to  determine  the  single  most  critical  link 
in  a  highway  system. 

The  following  sections  of  this  paper  briefly  present  Wollmer’s  method  for  finding 
the  most  vital  link  in  a  network  and  develop  an  improved  algorithm  for  solving  the  most 
vital  link  problem.  An  example  using  the  improved  method  is  included. 
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SECTION  II 

WOLLMER’S  ALGORITHM 


Wollmer’s  algorithm  for  finding  the  most  vital  link  in  a  network  follows  as 
a  consequence  of  the  following  theorem.  The  proof  of  this  theorem  is  given  in 
reference  (1). 

THEOREM  (Wollmer) 

Suppose  the  network  [N;A]  has  a  maximal  flow  value  of  v(F*),  while  the 
network  ^N;  jA-(x,y)jj  has  a  maximal  flow  value  of  v(F*jjy).  Then,  every  maximal 
flow  pattern  of  IN;A1  has  at  least  v(F*)  -  v(F*jjy)  units  of  flow  through  the  arc 
(x,y).  Moreover,  there  is  a  maximal  flow  pattern  of  (N;A]  which  has  exactly 
v(F*)  -  v(F*jjy)  units  of  flow  over  the  arc  (x,y). 

As  Wollmer  points  out,  “the  above  theorem  reduces  the  problem  to  one  of 
finding  that  link  whose  minimal  flow  among  all  maximal  flow  patterns  is  greatest.” 
Wollmer’s  iterative  procedure  for  finding  this  link  is  as  follows; 

Step  0: 

Find  a  maximal  flow  pattern,  F*,  in  the  network  IN;A)  and  let  f*(x,y)  be 
the  corresponding  flow  in  each  arc  (x,y)  e  A. 

Set  the  ‘least  flow’  of  each  arc  (x,y),  equal  to  P(x,y). 

Let  f*(p,q)  =  max  f*(x,y)  and  go  to  Step  1. 

(x.y)  f  A 

Step  1 : 

Solve  a  maximal  flow  problem  for  the  network  j^N;  jA  -  (p,q)  jj.  Let  the  cor¬ 
responding  maximal  flow  pattern  be  denoted  as  F*pq  and  go  to  Step  2. 

Step  2: 

(a)  Set  the  capacity  c(p,q)  of  arc  (p,q)  equal  to  v(F*)  -  ''(F*pq)  and  solve  a 
maximal  flow  problem  for  this  network  ^i.e.,  the  network  (N;A1  with 
c(p,q)  =  v(F*)  -  v(F*p^)j.  Call  the  corresponding  maximal  flow  pattern  F 
(Note  F’is  also  a  maximal  flow  pattern  of  IN;A] .) 

(b)  Next,  compare  the  least  flow  of  each  arc  (x,y)  with  f '(x,y)  and  if  f '(x,y)  < 
least  flow  of  arc  (x,y),  replace  the  least  flow  of  (x,y)  with  f ’(x,y).  Reset 
c(p,q)  to  its  original  value  and  go  to  Step  3. 
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Step  3: 

Let  U  =  max  j  least  flow  of  arc  (x,y)  . 

(x.y)t  A 

(a)  If  U  <  v(F*)- v(F*pj^)  terminate;  (p,q)  is  a  most  vital  link. 

(b)  If  U  >  v(F*>  -  v(F*p^^ );  find  an  arc  (p,q).  such  that  the  least  flow  of  (p,q) 
equals  U,  and  go  to  Step  I . 
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SIXTION  III 

AN  IMI’ROVI  I)  AlCORITHM 


Wollmcr's  .il^onlhit)  consulcrs  oach  arc  as  a  L'andidalo  I'oi  the  must  vital  link. 
However,  a  necessary  condition  is  employed  in  the  improved  algorithm  which  reduces 
the  niimiH’r  ol  arcs  that  must  Iv  considered  explicitly  as  candidates. 

IIIIOKI  M  I 

A  necessary  condition  for  an  arc  (a.h)  to  be  a  must  vital  link  is  that  fur  any 
maximal  flow  pattern  in  the  network  i  N:A| .  the  flow  in  arc  (a.b)  is  at  least  as  great  as 
the  How  over  every  arc  in  a  minimal  cut. 

rite  following  lemma  will  he  useful  in  the  proof  of  Theorem  I. 


LI  MMA  I 

If  t  .X.Xl  IS  a  minimal  cut  cimtaining  at  least  two  arcs  in  a  network  |  N:A|  and  if 
arc  tx.y » is  in  (X..Xl.  then  (X.Xi  -  tx.y)  is  a  minimal  cut  in  the  network  FnJa  -Ix.y)! 


PROOl  Ol  LI  MMA  I 

Suppose  that  ( Y.Y  I  is  a  minimal  cut  in  ^N;j  A  -  (x.y)  j  j  and  that  C(Y.Y)  < 
CtX.Xt  -  ctx.yf  Note  that  (Y.Y  tut  x.y)  has  to  be  a  cut  for  |N;A)  and  that  C(Y,Y)  + 
cl  x.y )  <  (’( X.X).  but  (X.X)  is  a  minimal  cut  of  |N;A|  and  (Y.Y)U(x.y)  is  a  cut  of  |N;A1 . 
Thus.  It  must  be  (rue  that  C'(  Y.Y )  *  ( (X.X)  -c(x.y)  and.  hence.  (X.X)  -  (x.y)  is  a  mini¬ 
mal  cut  in  I^N.  [a  -  (x.y  )|  j  . 

PK(M)| Ol  TIII  ORI  M  I 

Let  (X.X)  be  a  minimal  cut  in  |N:A|  and  note  that  by  hypothesis  arc  (a.b)  is  a 
most  vital  link  Assume  (hat  f*(a.b)<f*(p.q)  =  max  f*(x.y)  for  some  maximal  flow  pat- 

Is.VIclX.Xl 

tern  I  *  defined  in  |N.A|  l(  will  he  shown  that  this  assumption  leads  to  a  contradiction; 
(a)  H>  l  emma  I  v(  I  I  =  v(|-* )- f•(p.(J),  by  assumption  f*la, I'*  <  f*(p.(i); 

therefore.  *p,,  •  =  v(l  •)  -  (•(p.«|><  v(|  • )  -  f*(a.b). 

(bi  I  urilier.  (here  exists  a  (low  pattern  in  |^N.  |  A  -  a.b)  j  j  with  the  value 
vll  I  -  v(l  •  I  -  rta.h).  and  hence,  (he  maximal  How  value.  v(l'*j||,).  in 
fs  |A-(a.b)j  I .  must  satisfy  ^  v(l* )  -  (♦(a.b). 
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(c)  Conditions  (ii)  and  (b)  imply  that  >  v(F*pj^);  but  this  leads  to  a 

contradiction,  since,  by  assumption,  (a,b)  is  a  most  vital  link,  and  by  the 
definition  of  a  most  vital  link,  it  follows  that;  v(F*j,|,)  <  vlF*^^ ),  for  all 
<x.y)cA.  O.I.D. 

Theorem  I  guarantees  that  those  arcs  whose  flow  is  less  than  the  largest  How 
through  an  arc  in  a  minimal  cut  for  some  maximal  flow  pattern  in  |  N;A|  need  not  be 
considered  as  candidates  for  the  most  vital  link. 


SKTION  IV 
A  LABI  LIN(;  SCHI  Ml; 


A  labclin'.'  scheme  is  employcil  in  the  algorithm  for  finding  a  most  vital  link  from 
the  set  ot'caiulidafe  ares.  At  tlie  outset  of  the  labeling,  it  is  assunicd  that  there  is  a 
ma.ximal  How  pattern.  I  defined  for  the  network  |  N;A| .  Suppose  that  arc  (a.b)  is  an 
arc  from  the  candidate  set  (this  set  is  made  up  of  those  arcs  that  satisfy  the  necessary 
condition  of  Theorem  1 ).  Next,  an  attempt  is  made  to  label  from  node  a  to  node  b  without 
using  the  arc  (a.b).  The  labeling  (and  backtracking)  rules  to  Ih'  used  are  similar  to  those 
employed  by  Lord  and  h'ulkerson  in  their  algorithm  for  the  solution  of  the  maximal  flow 
problem.^'^*  The  labeling  and  How  changing  rules  are; 

Let  node  a  be  labeled  with  (-.  c  (a)  =  »»).  At  a  general  step,  suppose  the  node  x  is 
labeled  (/.i.  e  ( x ))  and  that  the  nodes  x  and  y  are  connected  by  some  arc.  Then,  node  y 
may  be  labeled  if  either  of  the  following  situations  occur: 

(a)  fty.x)  >  0;  then  node  y  is  labeled  (x~.  c  (y))  where  e(y)  =  min  (e(x).  f(y,x)); 
or 

(b)  f(x.y)  <  c(x.y );  then  node  y  is  labeled  (x'*’.  e  (y))  where  e(y)  =  min  (e  (x). 
c(x.y)  -  f(x.y )). 

The  labeling  process  is  continued  until  either  node  b  is  labeled  (breakthrough),  or 
node  b  is  unlabeled,  in  which  case  no  more  labeling  is  possible  (non-breakthrough).  If 
breakthrough  occurs,  node  b  must  have  a  label  of  the  form  (ip.  c  (b))  for  some  node  q. 
Likewise,  node  q  has  a  label  (r*.  e  (q ))  and  node  r  has  a  label  (p~.  e  (r)),  etc.  Thus  a 
series  of  nodes  (and  a  path  of  arcs)  starting  with  node  h  and  ending  with  node  a  is  defined. 

Let  f  =  min  |e  (b).  f(a.b)| .  Tor  each  arc  'x.y )  on  this  path,  change  the  flow  as 

follows. 

(a)  If  node  y  has  a  label  (x'*'.  6  (y)).  replace  ftx.y)  by  ftx.y)  +  f. 

(bl  If  node  y  has  a  label  (x”.  e  (y)).  replace  f(y.x)  by  f(y.x)  -  e. 

Tinally.  decrease  the  How  in  arc  (a.b)  by  t  units. 

As  l  ord  and  l  ulkerson  have  shown,  the  labeling  scheme  must  end  in  one  of  two 
mentioned  states,  breakthrough  is  achieved  or  breakthrough  is  not  achieved.  The  follow¬ 
ing  results  indicate  what  can  be  deduced  when  either  of  these  states  occur  at  the  end  of 
the  labeling  process. 


t) 
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THEOREM  2 


Let  F*  be  a  maximal  flow  pattern  in  the  network  (N;A] .  Let  (a,b)  be  an  arc  in 
(N;A1  and  use  the  labeling  scheme  to  label  from  node  a  to  node  b  without  using  arc 
(a,b).  If  breakthrough  occurs,  the  value  of  arc  (a,b)  is  less  than  or  equal  to  f*(a,b)  -  e. 

If  breakthrough  does  not  occur,  the  value  of  arc  (a.b)  is  equal  to  f*(a,b). 

PROOF 

If  breakthrough  occurs,  an  alternate  maximal  flow  pattern  is  found  by  making 
an  €  change  of  flow  in  the  path  established  by  the  labeling  and  decreasing  the  flow  in 
(a,b)by  c.  Then  v(F*,j,)>  v(F*)-  f*(a,b)  +  c,  and.  since  v(a,b)  =  v(F*)  -  v(F*,|,)  < 
v(F*)  -  ^v(F*)  -  f*(a,b)  +  cj ,  and  e  >  0,  it  follows  that  v(a,b)  <  f*(a,b)  -  c  <  f*(a,b). 

Assumethat  breakthrough  does  not  occur,  then  by  the  nature  of  the  labeling 
rules,  there  is  no  maximal  flow  pattern  in  (N;A1  with  less  than  f*(a,b)  units  of  flow  over 
arc  (a,b),  and,  hence,  by  Wollmer’s  Theorem.  v(a,b)  -  f*(a,b). 

O.E.D. 


Utilizing  these  results,  the  following  algorithm  can  be  employed  to  locate  a  most 
vital  link  in  a  network: 

Step  0: 

(a)  Find  a  maximal  flow  pattern,  F*,  in  the  network  |N;A)  and  let  (X,X)  be  a 
minimal  cut.  Let  U*  *  max  c(x,y)or,  alternately,  U*  =  max  f*(x,y). 

(>.y)<(X.X)  (».y)c(X.X). 

(b)  Note  those  arcs  (x,y)  c  A  for  which  f*(x,y)  >  U*  and  store  these  arcs  in  a 
list:  these  arcs  form  the  candidate  set,  S.  For  each  arc  in  this  set.define  an 
upper  bound  as  U(x,y)  =  f*(x,y). 


Step  I : 

Let  U(a,b)  =  max  U(x.y)  and  set  ftx.yt  =  f*(x.y)  for  all  arcs  in  |N;A|. 


Step  2: 

Use  the  labeling  rules  to  label  from  node  a  to  node  b  without  using  arc  (a.b). 
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Step  J 


(;i)  ir  hreukthroi'gh  occuri.  use  the  backtracking  and  flow  changing  rules, 
replace  Kx.y )  by  the  resultant  flow  in  each  (x.y)  c  A.  and  if  fta.b)  >  0, 
repeat  Step 

(h)  Otherwise,  replace  U(a.b)  with  fta.b)  and  if  U(a.b)  >  max  U(x.y)  terminate; 

I  •  III  (».y)»  s 

arc  ta.b)  is  a  most  vital  link. 

Otherwise,  replace  f*(x.y)  with  ftx.y)  for  all  (x.y )  e  A  and  go  to  Step  I . 

By  the  use  of  this  algorithm,  a  maximal  How  pattern  is  always  maintained  and  if 
there  is  no  breakthrough  for  the  arc  (a.b)  and  U(a.b)  >  Utx.y  I  for  all  arcs  (x.y)  of  the 
candidate  set  S.  the  arc  (a.b)  can  be  declared  a  most  vital  link.  i.e..  non-breakthrough 
implies  that  v(a.b)  -  f*(a.b)  =  U(a.b)  and.  thus,  v(a.b)  >  U(x.y )  >  v(x,y)  for  (x.y)  e  S. 
and  by  Theorem  I  the  most  vital  link  of  iN;A|  is  an  element  of  S. 

Tlie  other  aspects  of  the  method  that  must  lye  considered  are  (I )  is  the  method 
finite,  and  (2)  is  there  a  way  to  locate  alternative  optimal  solutions?  The  finiteness  of 
the  procedure  follows,  since  there  can  be  only  a  finite  number  of  arcs  in  the  cand'date 
set.  and  there  is  a  finite  number  of  labelings  that  can  be  made  for  each  arc  of  the  set. 
Alternative  optimal  solutions  are  readily  identified  as  follows:  after  an  optimal  solution 
has  been  found  and  other  candidate  arcs  exist,  reapply  the  algorithm  to  any  arc  in  the 
candidate  set  whose  upper  bound  is  equal  to  the  value  of  the  most  vital  link  (a.b),  i.e., 
set  U*  =  U(a.b).  delete  arc  (a.b)  from  the  candidate  set.  and  reapply  the  algorithm.  If, 
upon  reapplication  of  the  algorithm,  the  most  vital  link  has  a  value  equal  to  U*.  then  an 
alternative  solution  has  been  found.  In  the  latter  case,  the  algorithm  is  reapplied  using  a 
further  reduced  candidate  set.  and  the  procedure  is  continued  until  all  alternative  optimal 
solutions  have  been  found. 


*  An  alternative  to  Step  .la  is  to  lest  if  ftx.y  l<  (Kx.y).  (x.y)c  S.  and  if  ftx.y )  >  U*.  to  replace  Utx.y) 
willi  ftx.y).  il  It  x.y  I  <  I'*  then  arc  (x.y  I  may  he  dropped  from  llte  set  S.  lire  computations  are  con¬ 
tinued  hy  either  repealinii  Step  2  if  arc  (a.b)  has  not  been  dropped  or  returning  to  Step  1  if  arc  (a.b) 
has  been  dropped 
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SECTION  V 
EXAMPLE 


Consider  the  flow  network  (N;A|  shown  in  Figure  I. 


Figure  1.  Maximal  Flow  in  Sample  Network  with  Source  Node  >  and  Sink  Node  t 


Associated  with  each  arc  in  the  network  is  an  ordered  pair  of  numbers.  The  first 
number  corresponds  to  the  capacity  of  the  arc  and  the  second  number  corresponds  to 
the  amount  of  flow  over  the  arc  in  some  maximal  flow  pattern. 

The  improved  algorithm  will  be  employed  to  find  all  most  vital  links  in  the  net¬ 
work  I N  ;A  I . 

Step  0: 

(m)  Figure  I  presents  a  maximal  flow  pattern  F*  in  the  network  with  21  units  of 
flow  from  node  s  to  node  t.  The  individual  arc  flows  arc: 


0 


I*(s.a)  -  :5  |••(a.c)  =  15 

(■♦(s.b)  -  :  f*(a.tl)=  10 

=  17  l♦(b.c)=  2 

f*(d.t)  =  10  l*(b.d)=  0 


Ilic  corrospundirti!  niinitual  cut  (which  is  uni(|ui' I  contains  the  arcs  (a.c). 
(b.c).  and  (d.t).  Hence. 

U*  =  max  r*(x.y)  =  f*(a.c)  =  15. 

(i.ytf  (X.X) 


(b) 


Step  I : 


Tlie  set  of  candidate  arcs  is: 

S=  |(s.al.  (c.t).  (a.c) 


Start  with  arc  (s.a).  since  U(s.a)  =  f*(s,a)  =  max  LKx.y). 

(».y)f  s 


Step  2: 

Labeling  from  node  s  to  node  a  without  using  arc  (s.a)  is  possible  over  the  path 
containing  the  arcs  (s.b).  (b.d).  and  (a.d).  The  value  of  e  =  8. 

Step 

(a)  Since  breakthrough  occurred,  the  e  units  of  flow  are  removed  from  arc  (s.a) 
and  added  to  each  arc  in  the  path  found  in  Step  2.  The  revised  maximal 
flow  pattern  is  shown  in  Figure  2. 


Figure  2.  Revised  Maximal  Flow  in  Sample  Network 
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Step  2: 

Starting  with  the  maximal  flow  pattern  of  Figure  2.  it  is  not  possible  to  label 
from  node  s  to  node  a  and  non-breakthrough  has  occurred. 

Step  3: 

(b)  Replace  U(s,a)  with  the  flow  f(s.a)  =  1 7.  Since  U(s,a)  >  max  U(x,y )  = 

LKc.t),  arc  (s,a)  is  a  most  vital  link. 

Since  LKc.t) «  LKs.a)  =  1 7,  arc  (c.t)  may  also  be  a  most  vital  link.  In  order  to  test 
this  hypothesis,  arc  (s,a)  is  dropped  from  the  set  S  and  the  algorithm  is  reapplied  starting 
with  the  maximal  flow  pattern  in  Figure  2. 

Employment  of  the  labeling  rules  from  node  c  to  node  t  results  in  non-breakthrough, 
and  arc  (c,t)  is  an  alternative  most  vital  link. 
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